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Abstract 

Block-to-block interface interpolation operators are constructed for 
several common high-order finite difference discretizations. In contrast 
to conventional interpolation operators, these new interpolation oper- 
ators maintain the strict stability, accuracy and conservation of the 
base scheme even when nonconforming grids or dissimilar operators 
are used in adjoining blocks. The stability properties of the new opera- 
tors are verified using eigenvalue analysis, and the accuracy properties 
are verified using numerical simulations of the Euler equations in two 
spatial dimensions. 
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1 Introduction 

Wave-propagation problems frequently require farfield boundaries to be po- 
sitioned many wavelengths away from the disturbance source. To efficiently 
simulate these problems requires numerical techniques capable of accurately 
propagating disturbances over long distances. It is well know that high-order 
finite difference methods (HOFDM) are ideally suited for problems of this 
type. (See the pioneering paper by Kreiss and Oliger [I3]). Not all high-order 
spatial operators are applicable, however. For example, schemes that are G- 
K-S stable [8], while being convergent to the true solution as Ax — > 0, may 
experience non-physical solution growth in time [3], thereby limiting their 
efficiency for long-time simulations. Thus, it is imperative to use HOFDMs 
that do not allow growth in time; a property termed "strict stability" [7] . 
Deriving strictly stable, accurate and conservative HOFDM is a significant 
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challenge and that has received considerable past attention. (For examples, 
see references [ISlESlEZlIlHlEElIiniEnii). 

A robust and well-proven HOFD methodology that ensures the strict sta- 
bility of time-dependent partial differential equations (PDEs) is the summation- 
by-parts simultaneous approximation term (SBP-SAT) method. The SBP- 
SAT method simply combines finite difference operators that satisfy a summation- 
by-parts (SBP) formula P^. with physical boundary conditions implemented 
using either the Simultaneous Approximation Term (SAT) method p], or the 
projection method [271 EHl HH]. Examples of the SBP-SAT approach can be 
found in references [2l[23l[25l[l8l[20l[22l[26l[I71[33l[3lE^ 

An added benefit of the SBP-SAT method is that it naturally extends to 
multi-block geometries while retaining the essential single-block properties: 
strict stability, accuracy and conservation [4J. Thus, problems involving com- 
plex domains or non-smooth geometries are easily amenable to the approach. 
References pTj [20l [35] report applications of the SBP-SAT HOFD methods 
to problems involving non-trivial geometries. 

Current multi-block SBP-SAT methods suffer from two significant im- 
pediments: 1) the collocation points in adjoining blocks must match along 
block interfaces (i.e., "conforming" grids must be used), and 2) identical SBP 
schemes must be used tangential to the block interfaces. These impedi- 
ments prevent SBP-SAT operators from use on adaptively refined noncon- 
forming multi-block grids and from use where hybrid-approaches involving 
two or more discretization techniques is desirable. 

Although neither impediment (i.e., conforming grids or identical ele- 
ments) significantly limits Finite- and Spectral-Elements methods, @ to date, 
mitigating their influence in the context of HOFDM has rarely been success- 
ful. A noteworthy exception appears in reference [53] , where a hybrid method 
is developed (based on the SBP-SAT technique) to merge a HOFDM to an 
unstructured finite-volume method. This was accomplished by alternating 
the finite volume scheme at the interface boundary. 

The goals of the present study are two-fold. The first is to develop a sys- 
tematic methodology for coupling arbitrary SBP methods between adjoining 
blocks. This task essentially requires identification of the conditions that 
adjoining interpolation operators must satisfy to guarantee interface strict 
stability, accuracy and conservation. 

The second task is to demonstrate the new methodology by developing 
interface interpolation operators that couple commonly used SBP operators. 

^Ironically, these impediments were assumptions needed to prove stability and conser- 
vation, in the original single domain SBP-SAT approach. 

^For example, see the nonconforming spectral multi-domain work of Kopriva and Kolias 
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To this end, new interpolation operators are derived to couple 2nd-:2nd-, 4tli- 
:4tli-, 6tli-:6tli-, and 8tli-:8tli-order SBP operators across a nonconforming 
interface containing a 2 : 1 grid compression, as well as to couple 2nd-:4th- 
and 4tli-:8tli-order SBP operators across a conforming interface. 

In Section [2] we introduce some definitions and discuss the SBP property. 
Sufficient conditions for the interface stability of a two-dimensional hyper- 
bolic system are introduced in Section [HJ The new approach utilizes the SAT 
technique in combination with newly developed SBP-preserving interpolation 
operators. The construction procedure for the SBP-preserving interpolation 
operators are presented in detail. In Section H] the accuracy and stability 
properties of the newly developed multi-block interface coupling are tested 
by performing an eigenvalue analysis and numerical simulations of an ana- 
lytic Euler vortex. In Section [3] conclusions are drawn. Sufficient conditions 
for interface stability of a 2D parabolic system followed by SBP-preserving 
interpolation operators are presented in the Appendix. 



2 Definitions 

Two- and three-dimensional schemes are constructed using tensor products 
of one-dimensional SBP finite-difference operators. Thus, we begin with a 
short description of ID SBP operators including some relevant definitions. 
(For more details see [T^ISlj and [I'S]). 

2.1 Summation-By-Parts 

Assume the existence of real- valued functions u,v such that u,v & L^[0, 1]. 
Define an inner product to be {u,v) = J^uvdx, and let the corresponding 
norm be = {u,u). With these definitions, consider the hyperbolic scalar 
equation Ut + = (excluding the boundary condition). Apply the energy 
method to Ut + Ux = 0; i.e., multiplication by u and integration by parts 
leads to ^ 

— II^IP = -{u,u,x) - {ux,u) = -u^lo , (1) 

where m^Iq = u^{x = 1) — u^{x = 0). 

Next, discretize the domain (0 < x < 1) using N-l-1 equidistant grid 
points, 

Xi=ih, i = 0,l...,N, h = . 

Denote the projection of v onto the discrete grid point and define 

the discrete solution vector to be v'^ = [vo,Vi, ■ ■ ■ ,vn-]. As with the contin- 
uous inner product, define the discrete inner product for real- valued vector 
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functions u,v E R^^"*" to be {u,v)h = Hv, where H = > 0. The 
corresponding discrete norm becomes ||f |||^ = Hv. 

Definition 2.1 A difference operator Di = H~^Q approximating d / d x is a 
first- derivative SBP operator if H = , x^Hx > 0, x ^ 0, and Q + = 
B = diag {-1,0..., 0,1). 

A semidiscretization of + = is Vt + Div = 0. 

A semi-discrete energy for the system can be derived in a manner analo- 
gous to that used in the continuous case. That is, multiplying Vt + Div = 
from the left by v^H and adding the transpose leads to, 

j^Ml = -{v, H-'Qv)h - {H~'Qv, v)h = -v^iQ + Q^)v = vl - vl . (2) 

Equation ([2]) is the discrete analog of ([1]). 

Three definitions are central to the present study. 

Definition 2.2 An explicit pth-order accurate finite difference scheme with 
minimal stencil width of a Cauchy problem, is called a pth-order accurate 
narrow stencil. 

Definition 2.3 A pth-order accurate narrow stencil SBP-operator is called 
a narrow- diagonal SBP operator if H is diagonal. 

Definition 2.4 Let the row-vectors x^f and x^c be the projections of the 
polynomials x^ onto equidistant 1-D grids corresponding to a fine and coarse 
grid respectively. We say that Ip2C one? Ic2f are pth-order accurate inter- 
polation operators if = Ip2cy^ f — x'^^ o,nd = Ic2F^c ~ ^ f vanish for 
k = 0..p — 1 in the interior and for k = . . . {p — l)/2 at the boundaries. 

Remark The boundary closure for a pth-order accurate narrow- diagonal 
SBP operator is of order p/2 (see [IH])- The convergence rate for narrow 
stencil approximations of fully hyperbolic problems (e.g., the Eulers equa- 
tions) drops to {p/2 -\- l)th-order. (See refs. P, EE] for more information on 
the accuracy of finite difference approximations). 

Remark We state here without proof that the accuracy of the pth-order 
narrow-diagonal SBP operator is preserved when using the pth-order accurate 
interpolation operators. Numerical experiments presented later herein, are 
consistent with this conjecture. 
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Figure 1: Domain 2-D 



2.2 Two-Dimensional Domains 

We begin by introducing the Kronecker product 



Co,g-l D 



Cp-lfiD ■■■ Cp^i^q-i D 



where C is a p x g matrix and D is an m x n matrix. Two useful rules for 
the Kronecker product are {A®B){C0D) = (AC) O (BD) and (A^B)^ = 
A^ ® 5^. 

Next, consider the 2-D rectangular domain Vi^^^ bounded by Z < x < 
f-i d < y < u with an (A^ + 1) x (M + l)-point equidistant grid defined as, 



0,1...,A^, 



yj=jhy^ j = 0, l...,M, 



h = — 
h = — 



The numerical approximation at grid point (xj, is a 1 x /c- vector denoted 
The tensor product derivations are more transparent if we redefine 



the component vector "vector of vectors". Specifically, define a 

discrete solution vector v'^ = v^, . . . , v^], where = [vp^, . . . , fp,M] 
is the solution vector at Xp along the y-direction, illustrated in Figure [U To 
distinguish whether a difference operator P is operating in the x— or the 
y-direction we use the notations and P^. The following two-dimensional 
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operators are frequently used, 



el =(eo®P), 

where D^'^ , and H^^y are the one-dimensional operators in the x- and y- 
direction, respectively, are identity matrices of appropriate sizes, and cq, 
and cn are one- dimensional "boundary" vectors defined by 

eo = [1,0,...,0]^, e;v = [0,...,0,l]^ . (4) 



Dy = [r ® Dv) 

el = {P ® eo) 
eL = {P (g) Cm) 



3 Analysis 

Our main focus is on multi-block interface coupling solving PDEs on non- 
conforming grids, combining different (here referring to the formal order of 
accuracy) finite difference SBP schemes. 

Consider the following 2-D hyperbolic problem (the extension to parabolic 
problems is shown in the Appendix): 

Vit + Au^ + Bviy = 0, X G 

V, + Av,. + = 0, ^enx^ ' 

where A and B are symmetric kxk matrices (fc = 4 for the compressible Euler 
equations). (The extension of the present study to 3-D is straightforward.) 
A corresponding semi-discrete approximation of ([5]) can be written 

Ut + A® D^M + 5 ® T>lu = SATj^ + SATl 
Vt + A® D^i; + 5 ® W^v = SATj^ + SAT^ , 

where SAT[ ^ denote the penalty terms at the outer boundaries. (These 
are tuned to obtain stability, see [211 ESI EH] for details.) The penalty terms 
handling the multi-block coupling can be written as 

SATIl = ® (H£)-ie^ {uN -Ik® If2cVo} 
SATlj, = J:R0CH.l)-'el{vo-Ik0lc2FUN} ' 

(However, the penalty terms (171) is not the most general representation, as 
will be shown later in this section.) The continuity condition u — v = along 
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the interface (at Xj) is approximated by utv — /a;(8>/f2C^o and VQ — Ik®Ic2F'UN 
in the above penalty-terms. and are unknown k x k matrices to be 
determined below by stability. We use the subscripts L, R to distinguish the 
operators in the left and right domain (which can be completely unrelated, 
i.e., we can use different numerical schemes in the two domains), with the 
only assumption that they fulfill the SBP property. If2C and Ic2F are 1-D 
interpolation operators constructed such that Ik ® If2C'Vo and Ik ® Ic2fUn 
have the same dimension as mat and vq respectively, un and vq are km x 1 vec- 
tors with components aligned with the block interface. For future reference 
we introduce the following norms: 



Ml = {Ik® HI® Hi) , MR = {h®H^R®Hl) 
and the following terms: 



jrpi 



-M^(A ® HI)UN 

+2«^(El ® Hl)uN 



-u 



R 



HIIf2c)vo 



^C2F 



+v^{A ® Hl)vo 
+2vI{J:r®HI)vo 

(E/j ® HI^Ic2f)un 
-Voi^L ® If2cHI)un 



(8) 



(9) 



X 



HI 

-H^Ic2F 



-H\If2C 

HI 



T 


UN 


,W = 









We introduce the following important relations. 



(10) 



HrIc2F 



Hii^L ~ If2cIc2f) > , Hl^{Il^ — Ic2fIf2c) > 0, 



(11) 

(12) 



where ^ are identity matrices of appropriate sizes. The following definition 
is central to the present study. 

Definition 3.1 A pth-order accurate interpolation operator that fulfills fTl\} 
is called a pth-order accurate SBP-preserving interpolation operator. 

The following lemma is central to the present study: 

Lemma 3.2 The scheme ^ with the penalty terms given by ([^ is stable 
if 2I1f = A , 2E/J = —A and 177]) hold, assuming that the boundary terms 
SAT[ R are correctly implemented. 
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Proof Apply the energy method by multiplying the first and second equation 
in by vF Ml and v^Mr (defined by ([H])) respectively, yielding, 

j^{\\u\\m, + \\v\\m,) = BT' + IT\ 

where BT^ corresponds to the outer boundary terms (which we assume are 
correctly implemented, i.e., bounded) and IT^ corresponds to the multi-block 
interface coupling given by (EI). 

A stable coupling is obtained if IT^ is non-positive. A non-dissipative 
and stable coupling is achieved by choosing 

2Sl = A , 2Sk = —A , H^Ic2F = -^f2C-^l ) 
which makes IT^ = 0. □ 

Remark Note the critical relationship that equation ffTTl) imposes on the 
two interface interpolation operators. Indeed, the two conditions 21]^ = 
A , = —A and equation ffTTl) provide sufficient conditions for interface 
stability. 

Remark Another important interface property is conservation, especially 
for problems with discontinous solutions. Indeed, the Lax-Wendroff theo- 
rem relates discrete conservation to convergence to the weak solution of the 
differential equation. It is shown in [U [23] for conforming grids that con- 
servation provides a necessary condition for interface stability. Likewise, the 
conditions for conservation are fulfilled by imposing the necessary conditions 
for stability in Lemma [3.21 (or Lemma [3. 3p . 

For non-linear problems (such as the compressible Euler equations) charac- 
teristic boundary conditions are often imposed (since they introduce damping 
[20] at the boundaries) which can also be utihzed at block interfaces (see for 
example [201 [SS]). Since A is a symmetric matrix, it can be diagonalized by 
R^AR = A, where R is a matrix consisting of the eigenvectors of A. We 
introduce the notation A = i?|A|i?^, and 2A'^ = A ± A. (The matrix B is 
also symmetric and can be rotated to diagonal form with a similar transfor- 
mation.) To couple characteristic variables using the SAT technique [20l [35] 
at the block- interface we choose S^, = A~ , S/j = —A^. 

The following lemma states the stability conditions for a characteristic 
multi-block coupling: 

Lemma 3.3 The scheme ^ with the penalty terms given by ^ is stable 
if Hl = A~ , = —A^, ( flT]) and ( fl^) hold, assuming that the boundary 
terms SAT[ ^ are correctly implemented. 
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Proof Apply the energy method by multiplying the first and second equation 
in by vF Ml and v^Mr (defined by ([H])) respectively, yielding, 



d 
di 



(Mm, + \MmJ = BT' + IT' , 



where BT' corresponds to the outer boundary terms (which we assume are 
correctly implemented, i.e., bounded) and IT' corresponds to the multi-block 
interface coupling given by ©. 

A stable coupling is obtained if IT' is non-positive. If = A~ , S^j = 
-A+and H^Ic2F = Ihc^l hold, IT' = -A® w^Xw where and X 
are given by ( ITOl) . (By assumption H^j^If2C = {HrIc2f)^ which makes X 
symmetric.) Hence, stability follows if X is positive semi-definite, since A 
is a symmetric and positive semi-definite matrix. Sylvester's theorem states 
that X > if the following conditions hold: 



By definition Hi > 0, > 0. Hence stability follows if holds. □ 

The extra conditions in fll2l) . besides the necessary conditions in fllip 
introduce a significant restriction in the construction of the interpolation 
operators (as will be shown in the coming section). The energy estimate for 
the characteristic coupling (as described in Lemma 13. 3p gave us an idea how 
to remove the extra conditions ffT^ and yet introduce a damping mechanism 
by the penalty coupling. Consider the following energy estimate 



where Q is symmetric and positive definite and a is a positive integer {w'" and 
X are given by (llOp ). For Q = A and a = 1 we recover the energy estimate 
from the characteristic coupling, presented in Lemma 13.31 By choosing a an 
even integer, stability follows if X is symmetric (which is true if condition 
(ITT]) holds), regardless of the extra conditions in ( fT2l) . We introduce the 



HI > 0, HI- Hlh^c{Hl)-^HlIc2F > 
iJ|>0, Hl-HlIc2F{Hir'HllF2C>0 



di 



{\Hm, + Mmu) = BT'-n® w^X^w , 



(13) 
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following penalties: 

SATIl = + Y ® (HD"^e^ {un -h® If2cVo} 

- ® (H^)-^lF2cH^e^ {Ic2fUn - Vq} 

A ' 
SATI^ = ~ ^ ® (HR)"^e^ {vo - h ® Ic2fUn} 

-n® (H^)-iH^e^ {vo -h® Ic2fUn} 
-VL® (H^)"^Ic2FH^eS {If2cVq - un} 

The first two penalty terms in (JT^ (,i-e., the first in SATjj^ and the first 
in SATj^ ) are identical to the non-dissipative penalty coupling defined in 
Lemma [2121 The following lemma resolves the issue with the extra conditions 
in (II2D, 

Lemma 3.4 The scheme ^ with the penalty terms given by is stable 
if Q is symmetric and positive definite and 177]) holds, assuming that the 
boundary terms SAT[ ^ are correctly implemented. 

Proof Apply the energy method by multiplying the first and second equation 
in (El) by uF Ml and Mr (defined by (JED) respectively, leading to ( |T3l) . with 
a = 2. BT^ corresponds to the outer boundary terms (which we assume are 
correctly implemented, i.e., bounded). Stability follows if VL is symmetric 
and positive definite. □ 

Lemma 13.41 is validated numerically in Appendix [H where we also show that 
condition ( ITTl) is necessary to guarantee a stable characteristic coupling. 

Remark For non-linear phenomena the addition of artificial dissipation (AD) 
is often required for stability. The addition of SBP preserving AD [21] does 
not alter the stability conditions in Lemma 13.21 Lemma 13.31 or Lemma 13.41 

3.1 Construction of Interpolation Operators 

This section briefly describes the construction of the SBP-preserving inter- 
polation operators If2C and Ic2F- In the present study we will consider a 
two to one ratio between the coarse and the fine grids along the common 
interface. (The analysis is not restricted to a two to one ratio, it merely 



10 



simplifies the construction.) The structure (here showing the upper part) of 
the interpolation operator Ip2C is given by: 



ai,r \ 



aq,i 

a2s 0'2s-l Q-2S-2 

a2s 



a2s 



0'2s~-2 0-28-1 023 

OQ ■ ■ ■ 023-2 0-23-1 023 



(The lower part of If2C is obtained by a permutation of both rows and 
columns.) The other interpolation operator is given by /c2F = {Hf^)~^{H^lF2c) 
For example, in the fourth-order case (see the Appendix) we have used 
g = 3, r = 11 and s = 2 , which leads to the following form: 



IC2F 



(K 



\ 



bii,s 

2a3 2ai 2ai 2a^ 
2a4 2a2 2ao 2a2 
2a3 2ai 
2a4 2a2 



\ 



2a4 

2ai 2a3 

2ao 2a2 2a4 



(The lower part of Ic2F is obtained by a permutation of both rows and 
columns.) Notice that the unknowns in Ic2F are completely defined by the 
unknowns in If2C and the known coefficients in and if^. 

Remark A pt/i-order accurate interpolation operator will preserve the order 
of the finite difference scheme (E]) utilizing narrow-diagonal SBP operators, 
i.e., resulting in (p/2+l)t/i order accuracy. SBP discretizations not restricted 
(see for example [T8l [3T] ) to a (p/2)t/i-order accurate boundary stencil, require 
a more accurate boundary closure of the interpolation operators to maintain 
the convergence rate. We omit a detailed study (see [36] for more infor- 
mation on the accuracy of finite difference approximations) of the accuracy 
requirement. 

The following procedure describes the construction of the pth-order ac- 
curate interpolation operators, given the y-norms in the left {Hf^) and right 
(if^) domains: 
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• Build If2C with the form given above. 



. Build Ic2F = {Hl)-Hl^cHl given Hi and Hi. 
• Solve for accuracy (see Definition \2A\ . 

(In step 2 above we have assumed that the left domain corresponds to the 
coarse grid, if not simply replace the positions between Hi and Hi if we 
have the opposite scenario.) If there are free parameters left after imposing 
the accuracy conditions they are tuned to obtain an optimal L2 error. For 
the pt/i-order accurate case this means that = (ej)^ ■ + (e^)"^ ■ (see 
Definition 12. 4p is optimized for k = (p— l)/2 + l... (resulting in linear 
equations to be solved) until there are no free parameters left. 

With these assumptions there exist by construction SBP-preserving inter- 
polation operators for the second, fourth, sixth and eighth order cases (based 
on narrow-diagonal SBP schemes), by using the symbolic mathematics soft- 
ware Maple, (see the Appendix). We have also constructed interpolation 
operators coupling fourth to second and eighth to fourth-order accurate SBP 
discretizations. 

Remark This is the first time (to our knowledge) that stability for multi- 
block coupling of finite difference schemes with different order of accuracy in 
more than 1-D have been shown. 

For the second-order and fourth-order cases the constructed interpolation 
operators fulfill the extra conditions in ( fT2l) . besides the necessary conditions 
in ( ITTl) . It is an open question if the conditions in (fT2l) can be met for the 
sixth and eighth-order cases. 



4 Computations 

We perform some numerical tests to verify the accuracy and stability prop- 
erties of the SBP-preserving interpolation operators. 



4.1 Numerical Validation of Lemma 13^2 



To test the stability properties numerically we can either perform : 1) a 
long-time integration, or 2) an eigenvalue analysis. We chose the latter in 
the present study. Consider with the following setup: 



A = B 



(0) 
(0) 



u 



,(0)' 



u 



= : at the outer boundaries 
(1) = ■ at the block-interface 



'1 " 




-^(0)- 




-1 


, u = 




V = 



(15) 
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Figure 2: The spectrum of the discrete approximation to (|T5l) multiphed by 
hii for the fourth-order discretization, with stable SBP-preserving interpola- 
tion (left) and unstable non-SBP interpolation (right). 15^ unknowns in the 
left (coarse) block. (Notice the different scales.) 

Despite it's simplicity, this model problem is very challenging to solve nu- 
merically (without introduction of artificial damping) since it is energy- 
conserving, i.e., the eigenvalues are purely imaginary (see [I6] for details)). 
This makes it an ideal benchmark problem testing for stability of a certain 
numerical scheme, including boundaries. 

The eigenvalues of the numerical approximation of the continuous test- 
problem (given by (JSj) and (ITSll ) employing SBP-preserving interpolation op- 
erators fulfilling the conditions in (ITTl) for the fourth (using 15^ unknowns 
in the left (coarse mesh) region) and the sixth (using 21^ unknowns in the 
left (coarse mesh) region) order accurate cases are shown in Figures [2] and 
[3] respectively. We have multiplied the eigenvalues with the grid-size {h^i) 
of the fine-mesh domain. (The interpolation operators are presented in the 
Appendix.) To further indicate the importance of Lemma [3l2l we also present 
the corresponding results employing interpolation operators that do not ful- 
fill the conditions in fllip (and hence a violation to Lemma [3.2p . referred to 
as non-SBP interpolation. 

The non-SBP interpolation operators as compared to the corresponding 
SBP-preserving interpolation operator fulfilling the conditions in (ITTl) differ 
only at the boundaries. The boundary closures of the non-SBP interpola- 
tion operators (presented in the Appendix) are defined by the following two 
requirement: 1) to have the same formal accuracy at all grid-points, and 2) 
to use a minimal stencil width. The non-SBP interpolation operators clearly 
introduce an instability due to eigenvalues to the right of the imaginary axis. 
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Figure 3: The spectrum of the discrete approximation to (1151) multiphed by 
for the sixth-order discretization, with stable SBP interpolation (left) and 
unstable non-SBP interpolation (right). 21^ unknowns in the left (coarse) 
block. (Notice the different scales.) 



4.2 Vortex in Free Space 

To test the accuracy of the present method a 2-D Euler-vortex (an analytic 
solution [20] to the compressible Euler equations) is run across a multi-block 
interface. The problem consists of 2 blocks (5x5 unit area) having non- 
matching gridlines. The blocks are patched together according to ([6]). We 
set the Mach number Ma = 0.3. The vortex is initiated at x=4 (one unit 
to the left of the interface at x=5) with a 10 degree angle of the background 
free-stream and then propagated to t = 1, see Figure HI 
The convergence rate is calculated as 

where w is the analytic solution and ^('^i) the corresponding numerical so- 
lution with grid size hi. \\w — w^'^'^'^Wh is the discrete I2 norm of the error. 
Since the solution consists of 2 blocks we add together the two different I2- 
errors. The standard explicit 4th-order Runge-Kutta method is used for time 
integration. 

The convergence results going from coarse to fine grid using second, 
fourth, sixth and eighth-order accurate schemes are presented in Tables [TH 
(The number M is the number of unknowns in the y-direction in the coarse 
grid.) 

The potential gain of a hybrid approach (here referring to the multi-block 
coupling of a higher order scheme on a coarse-grid domain to a lower order 
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Figure 4: Pressure contours for the vortex at t=0 (initial data) and at t = 1, 
computed with the sixth-order method using 41 x 41 grid points in the coarse 
(left domain). 
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Table 1: Z2-error and convergence rate q. Second order case. Coarse to fine 
grid. Ma = 0.3. Vortex centered at interface a,t t — 1. 
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101 


-5.53 


3.07 


-5.12 


2.99 


-4.79 


3.13 


-4.70 


2.99 


151 


-6.07 


3.05 


-5.65 


2.99 


-5.33 


3.08 


-5.23 


3.01 



Table 2: /2-error and convergence rate q. Fourth order case. Coarse to fine 
grid. Ma = 0.3. Vortex centered at interface at t = 1. 



M 


log{p) 


q[p) 


logihiu)) 


q(M) 


logikiv)) 


q(v} 


logihip)) 


q{p) 


51 


-5.13 




-4.38 




-4.04 




-4.25 




101 


-6.27 


3.79 


-5.66 


4.25 


-5.30 


4.19 


-5.53 


4.24 


151 


-7.02 


4.24 


-6.44 


4.44 


-6.09 


4.48 


-6.29 


4.33 



Tabic 3: /2-error and convergence rate q. Sixth order case. Coarse to fine 
grid. Ma = 0.3. Vortex centered at interface aXt — l. 



M 


log{p) 


q[p) 


logikiu)) 


q(u) 


logihiv)) 


q[v) 


logikip)) 


qip) 


51 


-4.72 




-4.44 




-3.93 




-4.00 




101 


-6.07 


4.48 


-5.84 


4.66 


-5.39 


4.83 


-5.30 


4.30 


151 


-6.93 


4.87 


-6.67 


4.75 


-6.27 


5.04 


-6.13 


4.74 



Table 4: Z2-error and convergence rate q. Eighth order case. Coarse to fine 
grid. Ma = 0.3. Vortex centered at interface a,t t — 1. 
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M 


log{p) q^P^ log{h{u)) logMv)) q^^^ logMv)) 


qip) 


51 

101 

151 


-4.53 - -4.12 - -3.75 - -3.72 
-5.28 2.51 -4.81 2.30 -4.47 2.39 -4.45 
-5.68 2.23 -5.18 2.11 -4.85 2.13 -4.84 


2.45 
2.21 


Table 5 
to nne , 


: /2-error and convergence rate q. Fourth to second-order case, 
grid. Ma = 0.3. Vortex centered at interface at t = 1. 


Coarse 


M 


log{p) q^P^ log{h{u)) logiUv)) q^^^ logMp)) 


q(p) 


51 

101 

151 


-4.74 - -4.43 - -3.93 - -4.01 
-6.07 4.41 -5.82 4.60 -5.38 4.81 -5.29 
-6.89 4.69 -6.63 4.63 -6.25 4.96 -6.11 


4.25 
4.65 


Table 6 
to fine 


: /2-error and convergence rate q. Eighth to Fourth order case, 
grid. Ma = 0.3. Vortex centered at interface at t = 1. 


Coarse 



scheme on a fine-grid domain) is found by comparing the results (which 
are almost identical) in Tables [6] and HI In Table [6] we couple an eighth- 
order accurate SBP scheme in the left (coarse) domain to a fourth-order 
accurate SBP scheme in the right (fine) domain. In Table H] the eighth-order 
stencil is used in both domains. In Table Owe couple a fourth-order accurate 
SBP scheme in the left (coarse) domain to a second-order accurate SBP 
scheme in the right (fine) domain. A comparison of Tables [5] and [2] indicates 
the benefit of employing higher (than second) order accurate methods for 
wave dominated problems (the error of the second order schemes clearly 
dominates) . 

5 Conclusions 

Our approach have been to use SBP operators, the SAT technique and SBP- 
preserving interpolation operators to enforce the interface conditions in a 
stable and accurate way for general hyperbolic (and parabolic) problems, 
defined on nonconforming grids. The main objective was to construct inter- 
polation operators that combined the following desirable properties: 

• Stability by construction without interfering with the existing SBP 
schemes. 

• Maintaining the overall convergence rate. 

• Maintaining simplicity of the numerical scheme. 
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To achieve the three properties above, we have constructed interpolation 
operators with the following requirements: i) They preserve the summation 
by parts rule based on the underlying SBP scheme, ii) They have the same 
order of accuracy as the underlying SBP scheme, iii) They are of minimal 
width in the interior. 

Numerical computations for the 2-D compressible Euler equations cor- 
roborate the stability and accuracy properties and also show that a careful 
boundary treatment is necessary to guarantee stability of the interface cou- 
pling. 
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APPENDIX 



I Numerical Validation of Lemma 
[3:41 



In this section we verify numerically the necessity of the extra conditions 
(fT2l) that turned up in Lemma [3731 We will also verify the modified penalties 
(|T4l) that removed the necessity of conditions ( JT2l) . in order to obtain an en- 
ergy estimate, yet introducing a dissipative coupling (referring to the energy 
estimate ffT^ ). 

Consider ([U]) with the setup given by ffT^ . The scheme with the 
penalty terms given by ^ was proven stable if S^, = A~ , S/j = —A"*", 
assuming that both the conditions in ( ITT]) and ( |T2i) hold. However, for the 
sixth-order case condition ( fT2i) could not be met. This led us to consider 
the scheme (E]) with the penalty terms instead given by (1141) . where stability 
follows if Q is symmetric and positive definite. In the present study we use 
Q = A. The two different schemes (that differ only in how the SAT coupling 
is done) will be referred to as the characteristic coupling and the quadratic 
coupling. 

The eigenvalues of the two different numerical approximations (the char- 
acteristic coupling and the quadratic coupling) employing the SBP-preserving 
interpolation operators fulfilling the conditions in ffTTl) for the sixth- (using 
21^ unknowns in the left (coarse mesh) region) order accurate cases are shown 
in Figure [5] (compare with the corresponding non-dissipative coupling, pre- 
sented in Figure [3]) . We have multiplied the eigenvalues with the grid-size 
(hfi) of the fine-mesh domain. A closer look at the eigenvalues for the charac- 
teristic coupling reveals some eigenvalues with a positive real part, the largest 
of order 10~^. This is a strong indication that when employing a characteris- 
tic coupling the conditions in ( JT2l) is necesary. (The interpolation operators 
for the second and fourth order accurate cases do obey the extra conditions 
in ffT^ . We have verified, although not presented here, that the second- 
and fourth-order accurate cases does not have eigenvalues with positive real 
part for the characteristic couphng). The quadratic coupling (not required 
to obey ( |T2|) ) on the other hand results in eigenvalues having non-positive 
real parts, in agreement with Lemma [3. 4[ 



3.3 and Lemma 
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-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 -0.12 -0.1 -0.08 -O.OB -0.04 -0,02 0.02 



Figure 5: The spectrum of the discrete approximation to (fTSl) multiphed by 
hii for the sixth-order discretization, with the characteristic couphng (left) 
and quadratic couphng (right). 21^ unknowns in the left (coarse) block.) 

II Parabolic problems 

In this section we extend the analysis to include parabolic problems (such as 
Navier-Stokes equations) of the form 



vit + Avix + Buy = Cii 

^xy ^yx 



'd,ui 
i, r 



Vt + + BVy = Cu^xx + Ci2^^xy + C2lVy^ + C 22^ yy ■, X G fij „ 



(17) 



where A and B are symmetric k x k matrices. The continuity conditions at 
the interface are given by: u = v,Ux = Vx,Uy = Vy. Parabolicity requires 
that, 

x^{C + C^)x > , (18) 



where 



C 



Cn Ci2 
C21 C22 



A corresponding semi-discrete approximation of ffT7|) can be written 
ut + A® T>lu + B® T>lu = Pl + SATIl + SAT^i + SAT^ 



Vt + A® T>^v + B® D^t; = Pr + SAT' + SATj^ + SAT^ 



R 1 



(19) 



where 

Pl = Cn ® D£D£m + C12 ® D£D[m + C21 ® D^D£m + C22 ® D^D[m 
Pr = Cn ® ^n^RV + Cu ® D^D^i; + C21 ® D^D^w + C22 ® D^D^w 
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Cn^l ® (H^)"^e^ {(4 ® T>Iu)n - 4 ® /f2c;(4 ® D^t;)o} 

SiC2i ® (H^)-i(D[)^e^ {uN -Ik® If2cVo} 

CsiSi ® (H^)-ie^ {(4 ® D^n)^ - 4 ® /^2c(4 ® D^t;)o} , 

5^T,^ =S^Cn ® (H^)-i(D^)^e^ {t;o - 4 ® Ic2fUn} 

Cn^l ® (H^)"^e^^ {(4 ® D^t;)o - 4 ® /c2f(4 ® ^Iu)n} 

E^C2i ® (H^)-^(D^)^eS {t;o - 4 ® /c2fM7v} 

® (H^)-'e^, {(4 ® D^t;)o - 4 ® /c2f(4 ® D[n)^} , 

SATj^ and ^AT/r are given by Eq. (CD, or Eq. dUD). ^AT^^ and SAT^f^ are 
the SAT penalty terms handling the viscous multi-block coupling. SAT^ ^ 
denote the viscous penalty terms at the outer boundaries. (These are tuned 
to obtain stability, see [211 [221 ES] for details.) 

Lemma II. 1 The scheme ( I7P|) is stable if 



— , Yj\ — / , S| — , — / , 

— ^ 1 ~ ^ ^ 1 ~ ^ ' ^R ~ ^ ^ 1 

and Lemma lSTB. Lemma lSTSi or Lemma hold, assuming that the boundary 
terms SAT^ ^ are correctly implemented. S^'^ are symmetric matrices. 

Proof Apply the energy method by multiplying the first and second equation 
in (|T^ by u^M^ and v^Mn (defined by ([HI)) respectively, yielding, 

Jti\H\M, + II^IIa/J = fiT^ + IT' + JT^ + (JT^)^ . 

BT^ corresponds to the outer boundary terms (which we assume are cor- 
rectly implemented, i.e., bounded). IT' is given by ([9]) and corresponds 
to the inviscid part. IT^ + {IT^y" correspond to the viscous part of the 
multi-block coupling given by: 

+M^Cn(4 + Si + Yl){Hl ® D^n)^ 

-ulC^^ ® P(Si ^Hl + Yl® l'E2FHl)ih ® D^t;)o 

-vlc^^ ® iy{j:% ® HI + Si ® il^cHlWk ® ^Iu)n 

+u%C2i{Ik + Si + YDiHl ® BIu)n 

-vlC,,{h - - ^l){Hl ® D^t;)o 
-VIC21 ® ® + Si ® I^2cHl)ih ® D[n),v • 
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A stable coupling is obtained if IT^ is non-positive and if JT^ vanish. By 
assumption Lemma 13.2^ Lemma 13.31 or Lemma 13.41 hold which means that 
IT^ is non-positive. By choosing: 




(in combination with fill I) ) IT^ = 0. □ 

Remark The symmetric matrices T? and can be chosen arbitrarily for 
stability. (Note that S^'^ are not powers of E, but different parameters). 
However, the specific choice will affect the eigenvalues of (fT9l) . We use T? = 
= Ik/ 2 since it leads to a compact spectrum (we have no proof that this 
particular choice is optimal). 



Ill Interpolation Operators 
III.l Second-Order case 

The discrete norm is given by if = hdiag{^, 1, . . . , 1, |). The interpolation 
operators are given by: 



'F2C 



11 


1 


1 






20 


2 


20 






1 


1 


11 


1 


1 


40 


4 


20 


4 


40 






1 


1 


11 






40 


4 


20 








1 


1 








40 


4 



1 


1 






4 


40 






11 


1 


1 




20 


4 


40 




1 


1 


11 


1 


40 


4 


20 


4 






1 


1 






20 


2 



_ J_ 

40 
11 
20 



22 



and 



Ic2F — 



11 


1 


10 


10 


1 


1 


2 


2 


1 


11 


20 


10 




1 




2 




1 




20 



20 



1 




2 




11 


1 


10 


20 


1 


1 


2 


2 


1 


11 


20 


10 



J_ 

20 



'20 



11 


1 


10 


20 


1 


1 


2 


2 


1 


11 


10 


10 



III. 2 Fourth- to Second-Order Case 

Fourth to second-order, with second-order accurate interpolation. The upper 
left corner of If2C is given by: 



= 

ai,2 = 

11.3 = 

01.4 = 

«1,5 = 
«1,6 = 
ffll,7 = 
ai,8 = 

ai,9 = 

11,10 

ai,ii 



198814379483797276626589 
292154366866608903386020 

33159311983113441321394 
73038591716652225846505 

1971623072622604595228923 
16828091531516672835034752 

95808836247685639307767 
824906447623366315442880 

1754()7()19514()572218fi04399 



10517.1'i72()71<)7 
237()727(i()<)(i:iN 



J18')fi720 
'1741.'-{7 



14023409(i()95<)7227:!(i2r,2«9(i() 

55326079761477589117471 
2804681921919445472505792 

120787724660864994449209 
14023409609597227362528960 

176943055421257964368859 
10617667207197920521896720 

_ 32831043671234696789193 
14023409609697227362528960 

_ 88481S4S47S606230181361 
4949438686740197892657280 



12.1 = 

02.2 = 

12.3 = 

02.4 = 

«2,5 = 
0.2,6 = 
0.2,7 = 
«2,8 = 
0.2,9 = 
12,0 = 
02,11 



1342363899420690372183 
101394760863706442939854 

87193543695192148191232 
253486877134263607349635 

5578077558511726076198789 
19467792163911445044451968 



3895790315853711272834919 
16223160136592870870376640 

2244981324207232772622413 
12167371)11)2444053 ir,27S248() 

2135247()(i(i988()8!)li9:if<:i59 
16223160136592870870376640 

2744378984383218872921257 
48669480409778612611129920 

177648545875941908529241 
16223160136592870870376640 



46672983468620710671107 
12167370102444653152782480 

3876549717482762134187 
3244632027318674174075328 

_ 704697674731763751811223 
97338960819557226222269840 



03.1 = 

03.2 = 

03.3 = 

03.4 = 

03.5 = 

03.6 = 

03.7 = 

03.8 = 

03.9 = 
03,10 
03,11 



6011306S72161601487699 
738978692662598990917680 

17239549725384800082266 
184744673165649747729395 

13200899856172905887766503 
70941954495609503128087680 

3693069541319881367608217 
11823659082601583854681280 

1607455065911582798562599 
8867744:!ll')511878<)ll)l()()li() 

351431(i:i9:i()r)278487727:i7r)l 
11823659082601583854681280 

8869555969232379886041937 
35470977247804751564043840 

652819712407377597489161 
11823659082601583854681280 

12904128303803864069SS69 
1773548862390237678202192 

_ 212102769S2319643393399 
11823659082601583854681280 

_ 300890344598183663971671 
70941964496609503128087680 



The interior of If2C is given by: 



292764282684548619564311 
4949438685740197892657280 

27348461476516948621801 
824906447623366315442880 

138809444447592059191981 
618679835717524736582160 



_ 233575073382358527482521 

824906447623366315442880 

_ 139882799512953873064261 
824906447623366315442880 
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The upper left corner of Ic2F is given by: 



&1,2 ■ 

61,3 

^2,1 

62,2 

&2,3 

63.1 ■ 

63.2 ■ 

63.3 ■ 

63.4 ■ 

64.1 ■ 

64.2 ■ 
^4,3 
64,4 
65,1 
&6,2 

65.3 ■ 

65.4 ■ 
&6,B ■ 



198814379483797276626589 
206226611905841578860720 



1342353899420590372183 
20622661190584157886072 



6()113()fi5721filfi()1487699 



20(i22(i(ill<)()." 
16579(i55<)<ll55(^ 



:4ir,7S86()720 
'21Jl)(i()G97 



51556(i5297(i4(i()394715180 

10899192961899018523904 
12889163244115098078795 

8619774862692400041133 
51556652976460394715180 

1971623072622604595228923 
23757305691552949884754944 



5578077558511726076198789 
7919101897184316628251648 

13200899856172905887766503 
39595509485921583141268240 

14345449851542882358651239 
118786528457764749423774720 

1628750216210655868232039 
19797754742960791570629120 

1298596771951237090944973 
2199750526995643507847680 

3693069541319881367608217 
6599251580986930623543040 

1340074612349330482468249 
19797754742960791570629120 

1754070195140572218604399 
14848316057220593677971840 

2244981324267232772622413 
4949438685740197892657280 

1607455065911582798562599 
4949438685740197892657280 

6801662777932010900407069 
14848316057220593677971840 

292764282684548619564311 
2474719342870098946328640 



66.1 ■ 

66.2 ■ 

66.3 ■ 

66.4 ■ 

66.5 ■ 

67.1 ■ 

67.2 ■ 

67.3 ■ 

67.4 ■ 
67.5 

67.6 ■ 

68.1 ■ 

68.2 ■ 

68.3 ■ 

68.4 ■ 
68,6 ■ 
63,6 ■ 

69.1 ■ 

69.2 ■ 

69.3 ■ 



237072760963890638174137 
" 19797754742960791570629120 



71174902232936303979453 
■ 2199750526995643507847680 



1171438797684261626767917 
2199760626996643507847680 

11445178695736567846643629 
19797764742960791570629120 

27348461476616948621801 
412463223811683167721440 

55326079761477589117471 
3959550948592158314125824 

274437S9843S321S872921257 
19797754742960791570629120 

8869555969232379886041937 
19797764742960791570629120 

6854257176134739780148789 
19797764742960791570629120 

138809444447692059191981 

309339917858762368291080 

292764282684548619564311 
2474719342870098946328640 



120787724660864994449209 
19797754742960791570629120 

17764854S87S941908529241 
6699261680986930523643040 

652819712407377697489161 
6599261580986930623543040 

11446178696736667846643629 
19797764742960791570629120 

233575073382358527482521 
412453223811683157721440 



27348461476516948621801 
412453223811683157721440 

176943055421257964368859 
14848316057220593677971840 

45572983468520710671107 
4949438685740197892657280 

129041283038038640695669 
989887737148039578531466 



69.4 = 

69.5 = 

69.6 = 

69.7 = 

610.1 ■ 

610.2 ■ 

610.3 ■ 
610,4 
610,5 ■ 
610,6 
610,7 ■ 

611.1 : 

611.2 ■ 

611.3 ■ 

611.4 ■■ 

611.5 ■ 

611.6 ■ 

611.7 ■ 

611.8 ■ 



6801662777932010900407069 
14848316057220593677971840 



139882799512953873064261 
412453223811683157721440 



138809444447592059191981 
309339!)178r,87(i2:i(i82')l()8l) 



2927642821)84r)48(il!)5(i4:ill 
2474719342870098946328640 



32831043571234596789193 
19797754742960791570629120 



3876549717482752134187 
1319850316197386104708608 



21210275952319643393399 
6599251580986930623643040 

1340074612349330482468249 
19797754742960791570629120 

233575073382358527482521 
412453223811683167721440 

233575073382358527482621 
412453223811683167721440 

27348461476516948621801 
412453223811683157721440 

1504186273085288913082967 
118786528457764749423774720 

704697674731763751811223 
39596509485921683141258240 

300890344598183663971671 

39595509485921583141258240 

14345449851542882358651239 
118786528457764749423774720 



138809444447592059191981 
309339917858762368291080 

139882799512953873064261 
412453223811683157721440 

138809444447592069191981 
309339917858762368291080 

292764282684548619564311 
2474719342870098946328640 



III. 3 Fourth-Order Case 

The discrete norm is given hy H — h diag{^, ||, ||, ||, 1, . . .). 
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The upper left corner of If2C is given by: 



11,1 = 

ai,2 = 
ai,3 = 
ai,4 = 

0,1,5 - 
11,6 = 

ai,7 = 

11.8 = 

11.9 = 
<Il,10 

ai,ii 



37625155150141581259 
74136705821905530032 

109494(100531368265009 
15754f)4')')87ir,4')2.1i:il8 

598!)f<48(j21i:S(j!)r>l!l<l6355 
2016518:i!)8:i5r,8:i04i68704 

259928664!)16563226325 
2520647997944788021088 

256587334010286818055 
5041295995889576042176 

113845943840605140771 
2520647997944788021088 

537477213186037277591 
10082591991779152084352 

89579548577713740579 
2520647997944788021088 

33361750807060821603 
5041295995889576042176 

_ 45876275221828488395 
2520647997944788021088 



891498675213598216803 
20166183983568304168704 



12.1 = 

12.2 = 

12.3 = 

12.4 = 

12.5 = 

12.6 = 

12.7 = 

12.8 = 

12.9 = 
12,10 
12,11 



9465638066209876131 
2187032821746213135944 

647731()271B84()14677 
185:i417(i4r, 

2378l)!)lj:S7:S 



47(i:iS25()8 
542318389255 



69985()5()2!)5878820350208 

2109162960890226252711 
8748131286984852543776 

776008005836532793785 
17496262573969705087552 

469836672403155B2089 
8748131286984852543776 

1087641244818460786021 
34992626147939410175104 

220860994180249341801 
8748131286984852543776 

77236836126760009613 
17496262673969705087562 

_ 96386094557526072417 
8748131286984852543776 



1793510683939834799913 
■ 69986060296878820360208 



13.1 = 

13.2 = 

13.3 = 

13.4 = 

13.5 = 

13.6 = 

13.7 = 

13.8 = 

13.9 = 
13,10 
13,11 



9465638066209876131 
3187878350341937791376 

54391101655540111975 
7')(i9(i9r>87585484447844 

I(i2929ir)4(i27.'i5 11894441 
51006053605471004662016 

2010507950879293630905 
6375756700683875582752 

6208709747109065353983 
12751513401367751165504 

1804441319301200448807 
6375756700683876682762 

68210841472168154533 
26603026802736602331008 

395393460093074051961 
6376756700683875582752 

180548419983997949901 
12761613401367751 165504 

_ 55143363449566679649 
6375756700683875582752 

_ 21221519586950580219 
1186187293160488480612 



The interior stencil is given by: 

_ 10513333512638366307 n = 

"4 1186187293150488480512 ^ 10513333512638366307 ^ 32 

^ _J_ 2 296546823287622120128 265006822749707021207 

3 32 "0 593093646575244240256 
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The upper left part of Ic2F is given by: 



61,1 

&1,2 ■ 
bl,3 ■ 

62.1 ■ 

62.2 ■ 
^2,3 
63,1 
63,2 
&3,3 
&3,4 

64.1 ■ 

64.2 ■ 

64.3 ■ 

64.4 ■ 
&5,1 ■ 
65,2 
^5,2 
65,2 
&5,2 



37625156150141581259 
37068352910952765016 

556802239188816243 
18534176455476382508 

556802239188816243 
37068352910952766016 

1855830517480818051 
4633544113869095627 

6477310271684014677 
9267088227738191254 

921883078907459525 
9267088227738191254 

598984862036951996356 
25503026802735502331008 

23789903736542318389266 
25503026802735502331008 

1629291546275511894441 
25503026802735502331008 

515153342119279949043 
25503026802735502331008 

5304666630950269925 
74136705821905630032 

43044142058984209239 
74136705821905530032 

41030774507740686345 
74136705821905530032 

_J_ 
16 

85529111336762272685 
2372374586300976961024 

258669336278844264695 
2372374586300976961024 

2069569915703021784661 
2372374586300976961024 



171717780706426649681 

2372374586300976961024 

10513333512638366307 
593093646575244240266 



06,1 
&6,2 
&6,2 
&6,2 
^>6,2 
&7,1 
&7,2 
67,3 
67,4 
&7,5 
&7,6 
&8,1 
&8,2 
&8,3 
^8,4 
^8,5 
^8,6 
69,1 
fig, 2 



37948647946868380257 
1186187293150488480512 

16661222413438617363 
1186187293160488480612 

601480439767066816269 
1186187293160488480512 

147 
266 

_ J_ 
16 

537477213186037277591 
14234247517805861766144 

362547081606153695007 
4744749172601963922048 

19403613824056061511 
4744749172601953922048 

12985334314735644039143 
14234247517805861766144 

10513333512638366307 
148273411643811060064 

10613333612638366307 
593093646575244240256 

29859849525904580193 
1186187293150488480512 

73620331393416447267 
1186187293160488480612 

131797820031024683987 
1186187293150488480512 



147 

266 

_9_ 
16 

_ J_ 
16 

11120583602353607201 

2372374586300976961024 

25745611708920003171 
237237468630097696 1024 



69,3 = 
''9,4 = 
69,5 = 
f>9,6 = 
69,7 = 
610,1 ■ 
610,2 
610,3 ■ 
610,4 

610.5 ■ 

610.6 ■ 

610.7 ■ 
611,1 ■ 
611,2 
611,3 ■ 
611,4 
611,6 ■ 

611.6 ■ 

611.7 ■ 

611.8 ■ 



60182806661332649967 
2372374686300976961024 

171717780706426649681 
2372:i745S(i;i()()<)76t)61024 

265()()(iy227497()71J21207 
296546823287622120128 

10513333512638366307 
148273411643811060064 

10513333512638366307 
593093646676244240266 

_ 46876276221828488396 
3558561879451466441636 



32128698185842024139 
1186187293160488480512 

18381121149866659883 
1186187293150488480512 

49 
768 

_9_ 

16 

_9_ 
16 

_ J_ 
16 

2971662250711994(15601 
9489498345203907844096 

597836894646611599971 
9489498345203907844096 

304175114079624983139 
9489498345203907844096 



171717780706426649681 
9489498346203907844096 

10513333512638366307 
14827341 164381 1060064 

265006822749707021207 
296546823287622120128 

10513333512638366307 
148273411643811060064 

10613333612638366307 
693093646676244240256 



III. 4 Fourth-Order Case, Non-SBP Interpolation 

The upper left part of Ic2F' 



( 1 



_5_ 
16 



15 
16 

1 



_5_ 
16 



16 



2a3 2ai 2ai 2a3 

2a4 2a2 2ao 2a2 2a4 

2a3 2ai 2ai 2a3 

2a4 2a2 2ao 2a2 2a4 
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The lower left part of If2C'- 



( 1 



O4 03 O2 Oq Ol 02 04 

O4 03 02 Ol OO ^^1 0^2 O3 O4 



V 



ao-4 is listed in previous subsection. 



■ / 



III.5 Sixth-Order Case 

Sixth order accurate case is given below. Notice that the stability condition 
for the upwind formulation could not be achieved. The discrete norm is given 

\yu TJ — h H'innf^^^ 12013 2711 5359 7877 43801 1 \ 
uy n tLUlUy\^^^QQ, gg^^ , 4320 ; 4320 ' 8640' 43200' -*-'■■ V" 
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The upper left corner of If2C is given by: 



= 

ai,2 = 

11.3 = 
ai,4 = 
11,6 = 

«1,6 = 
a 1.7 = 
ai.8 = 
0-1,9 = 
11,10 

ai,ii 
ai,i2 
ai,i3 
ai,i4 

<Il,16 

11,16 
11,17 
a2.i = 
0.2,2 = 
0-2,3 = 

02.4 = 
02,6 = 

02.6 = 

02.7 = 

02.8 = 



493203253466037198549455362364169948587300751 
1056237669662369450269932613630339432431112782 



3494144 

65307081693884397225173946017702348750445730 
627618834831184725134966306765169716215566391 



114866096841096884636391939697630554903763548 
527618834831184726134966306766169716216666391 



58160488726161616986950806006548607936924127 
527618834831184726134966306766169716216666391 

15525 



I(i059678106898596175:-12()(i!)8()()48(i80817748g46115 
1055237669662369450269932613530339432431112782 



436768 

28603010322006197876985672016642805153037163 
527618834831184726134966306765169716215666391 

858196 
6988288 

42764668664746224527531090194001857246708700 
527618834831184725134966306766169716216656391 



6988288 

10622334830569558938004175701898426132770338 
627618834831184726134966306765169716215666391 



3494144 

_ 1167889404766086591410267452043339821491788 
527618834831184725134966306766169716216656391 

10324674409544267975653283320299930288476692 
464377248367170642761107058624806491383799467 

1201479 



63083470065854914766431215082154659242790665 
464377248357170642761107058624806491383799467 

721527 
3076328 

75656801964097460751067400512605712346209266 
464377248367170642761107058624806491383799467 

838645 
6160666 

30937093119827657161778923229206427591240999 
46437724836170642761107068624806491383799467 

345 



02.9 = 

02.10 ' 
(12,11 ■ 
"2,12 ■ 

12.13 ■ 

12.14 ■ 

12.15 ■ 

12.16 ■ 

12.17 ■ 

13.1 = 

13.2 = 

13.3 = 

13.4 = 

13.5 = 

13.6 = 

13.7 = 

13.8 = 

13.9 = 

13.10 ■ 

13.11 ■ 

13.12 ■ 

13.13 ■ 

13.14 ■ 

13.15 ■ 

13.16 ■ 

13.17 ■ 



37294205185814962316045245318486095186507918 
' 464377248357170642761107068624806491383799467 



3628149322 i(i.r)8Sli444947422r)97(i!)0(i91li:-{9S553827 
928754496714341285522214117249612982767598934 

55219 



22445907670907889073586224542480749734259694 
464377248357170642761107058624806491383799467 



6160666 

6433667362017414116320173560603825973179166 
"464377248367170642761107058624806491383799467 

5409 



578944702382543295705133726021669910745894 
464377248357170642761107058624806491383799467 

5162337204772128987826641660149965144238346 
104797029908956098603626174638462632102012849 



1388032 

23386090748321137420688072323840400038565838 
104797029908956098603626174638462532102012849 



120607102321976279545796687898247496634070403 
419188119635824394414504698653860128408061396 



694016 

14476182846428389369733457155566337017487527 
209594059817912197207252349276926064204025698 

217175 
1388032 

6622670052963322216326122365241056731782031 
209594059817912197207262349276925064204025698 



1388032 

21321625226448423063859202901016977066304133 
209594059817912197207252349276925064204025698 

78363 
694016 

8305592718705970541713249560330845208737268 
104797029908966098603626174638462632102012849 



1892889085160894469318781649971017262255710 
104797029908956098603626174638462532102012849 



1388032 

192981667460847765236044675340556636916298 
104797029908966098603626174638462632102012849 



The interior of If2C is given by: 



03 



321100777805071156797079160300427016498 
38656226451108852306759931626138890483959 

3 
512 

192K(iO-i66K830i2K!)i0782i71!X-il802562098988 
38656226451108852306759931626138890483959 

25 
512 



02 

ai 



4816511667076067351956187404506405247470 
38656226451108852306759931626138890483959 

75 
256 

25812195338906006034876765214121809824039 
77312452902217704613519863252277780967918 
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III. 6 Sixth-Order Case, Non-SBP Interpolation 

The upper left part of Ic2F'- 



( 



1 

63 



315 



105 



63 



45 



256 


256 
1 


128 


128 


256 


256 






7 

256 


105 
256 


105 
128 

1 


35 
128 


21 

256 


3 
256 






2a5 


2a3 


2ai 


2ai 


2a3 


2a5 






2a6 


2 04 


2a2 


2ao 


2a2 


2 04 


2a6 








2a5 


2a3 


2ai 


2ai 


2a3 


2a5 






2a6 


2 04 


2a2 


2ao 


2a2 


204 



The lower left part of Ipic'- 
( 1 

1 

1 

Clg 05 04 03 (i2 Ctl Cto ''^l *^2 ^3 '^4 '^5 ^6 

Oe 05 04 03 02 Oi Oo Oi Oji 03 04 05 

ao-6 is listed in previous subsection. 

III. 7 Eighth-Order Case 

The upper left corner of Ific is given by: 
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= 


0.49670323986410765203 


0-1,2 = 


0.80670591743192512511 


ai,3 = 


-0.014476656476698073533 


ai,4 = 


-0.19497555250083395132 


ai,5 = 


0.16074268813765884450 


11,6 = 


0.22100911221277072755 


11,7 = 


-0.53042418939472250452 


11,8 = 


-0.86254139670456354517 


ai,9 = 


0.0064231825172866668299 


«1,10 = 


= 0.69727164011248914487 


= 


= 0.61825742343094006838 


ai,i2 = 


= 0.43307239695721032895 


11,13 = 


= -0.15848187861199173328 


11,14 = 


= -0.83836831742920925562 


ai,i5 = 


= -0.41767239062238727391 


11,16 = 


= 0.39252899650233765024 


ai,l7 = 


= 0.33094736964907844809 


ai,i8 = 


= -0.073592865087013788440 


11,19 = 


= -0.099669762780958210515 


11,20 = 


= 0.010748160731101219580 


11,21 = 


= 0.018056833808814782786 


11,22 = 


= -0.00077275234787125894193 


11,23 ~ 


= -0.0014911993994710636483 


12,1 = 


0.0025487859584304862664 


12,2 = 


0.47007080279424540800 


12,3 = 


0.093589176759289320118 


12,4 = 


0.33386224041716468423 


12,6 = 


-0.11418395501435496457 


12,6 = 


-0. 18998278818813064037 


12,7 = 


0.38484431284491986603 



12,8 = 


= 0.65828018436035862232 


12,9 = 


= 0.039704539050575728856 


12,10 


= -0.47252462545568042557 


12,11 


= -0.44892698918501097924 


12,12 


= -0.34402935194949605213 


12,13 


= 0.087284452472801643396 


12,14 


= 0.59539401290875351190 


12,15 


= 0.30484430526276345964 


12,16 


= -0.27096308216867871783 


12,17 


= -0.23108785344796321535 


12,18 


= 0.050608234918774219796 


12,19 


= 0.068801842319351676214 


12,20 


= -0.0073322707316320135247 


12,21 


= -0.012333488857215226757 


12,22 


= 0.00052275043402033040521 


12,23 


= 0.0010087644967132781708 


13,1 = 


= -0.022656973277393401539 


13,2 = 


= -0.93771184228725468159 


13,3 = 


= -0.052699680965389826267 


13,4 = 


= 1.9581430555012001670 


13,5 = 


= 1.6586148101136731602 


13,6 = 


= 2.7435837109255836264 


13,7 = 


= -3.0164708462284474624 


13,8 = 


= -5.8235016129647971089 


13,9 = 


= -1.0472767830023645070 


13,10 


= 3.2615209891555982371 


13,11 


= 3.5478595826728208891 


13,12 


= 3.0921640208240511054 


13,13 


= -0.34364793368678654581 


13,14 


= -4.5566547247355317663 


13,16 


= -2.4329078834671892590 



13.16 = 1.9882413967858906122 

03.17 = 1.7259601262271516763 

03. 18 = -0.36895458453625989435 

03.19 = -0.50448363278283993228 

13.20 = 0.052797763052066775846 

13.21 = 0.088972343374038343284 

13.22 = -0.0037175165926095403240 
a3_23 = -0.0071737841052106668688 

04.1 = 0.0021626748627943672418 

04.2 = 0.027636709246281031633 

04.3 = -0.00055259484658035070394 

04.4 = -0.033553649469391355855 

04.5 = -0.049244245327794891233 
a4,6 = 0.093489376222103426899 
a4,7 = 0.45734620580442859300 
a4,8 = 0.66579609152388478842 

04.9 = 0.24716623315221892947 

04.10 = -0.15893464065423852815 

04.11 = -0.25881084331023040874 

14.12 = -0.26865808253702445366 
a4,i3 = -0.018060020510041282529 
a4,i4 = 0.30841043055726240020 

04.15 = 0.17712461121229459766 

04.16 = -0.12549015561536110372 

04.17 = -0.11250298507032009025 

04.18 = 0.022964117700293735857 

14.19 = 0.031709446610808019222 

04.20 = -0.0032149051440245356510 

04.21 = -0.0054335286230388551025 

14.22 = 0.00022177994574852164638 

14.23 = 0.00042797426992744435493 
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The interior of If2C is given by: 



as = 


-0.0023556211403911721282 




ar = 


-0.001220703125 


as = 


ae = 


0.018844969123129377026 


a2 = 


05 = 


0.011962890625 


ai = 


04 = 


-0.065957391930952819589 


ao = 



-0.059814453125 



0.33510652017261795103 



III. 8 Eighth- to Fourth-Order Case 

The upper left corner of Ip2C is given by: 
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ai,i = 0.58483823651420605520 
ai,2 = 0.65474734710702549940 
ai,3 = -0.083569673821672843574 
oi,4 = 0.022425968859768418924 
ai,5 = 0.31477269554917488790 
ai,6 = -0.19106598818802213021 
aij = -0.69007595018747002513 
ai,g = -0.30500476635452414633 
ai,9 = 0.011380303429615219398 
ai,io = 0.39023505096198663290 
ai,ii = 0.6929.3443221288252711 
ai,i2 = 0.27930147563078413713 
ai_i3 = -0.17509858064671991258 
ai_i4 = -0.36446581507281642455 
oi,i5 = -0.53431810490188879575 
ai_i6 = -0.047645255801647499040 
ai,i7 = 0.48950384943343111452 
ai,i8 = 0.14968600350317559693 
ai^ig = -0.19713593239026365497 
ai,20 = -0.051374514632076414476 
ai,2i = 0.048757624195804968082 
ai,22 = 0.0066939070511429668639 
ai,23 = -0.0055223124518961777586 

02.1 = 0.011443922859934377398 

12.2 = 0.37277008707234857179 
a2,3 = 0.35282527696092662640 
02,4 = 0.096066986437174930272 
02,6 = -0.22996810860646674853 

02.6 = 0.12123432064337482341 

02.7 = 0.49992261650980671917 



02.8 = 0.25040596791141593267 

02.9 = 0.039631571628602904441 



02,10 


= -0.25355792154009083332 


02,11 


= -0.50801440584689113247 


02,12 


= -0.21961026714024651967 


02,13 


= 0.095341708111114624790 


02,14 


= 0.24893964119173768791 


02,15 


= 0.39235409374512417332 


02,16 


= 0.042074861655253219822 


02,17 


= -0.34486664831340391286 


02,18 


= -0.10508512775619684648 


02,19 


= 0.13702066788245211707 


02,20 


= 0.035285372237949024781 


02,21 


= -0.033422056958152836350 


02,22 


= -0.0045282849361975826191 


02,23 


= 0.0037357262504306790489 


03,1 = 


= -0.095787414722590617699 


03,2 = 


= -0.68323289165328697626 


03,3 = 


= -0.43569930844225202810 


03,4 = 


= 1.8027331550089460412 


03,5 = 


= 3.7567210660509388125 


03,6 = 


= -0.040511780.395060609101 


03,7 = 


= -3.9354784185565657354 


03,8 = 


= -2.5330463966461275528 


03,9 = 


= -1.0862631313184671074 


03,10 


= 1.5476108406776833914 


03,11 


= 4.0781324522987606857 


03,12 


= 1.9733308525923243475 


03,13 


= -0.37565017560062195914 


03,14 


= -1.7659882182232900992 


03,16 


= -3.1675028227076565595 



03.16 = -0.41961367860804052329 

03.17 = 2.6138186452177357535 

03.18 = 0.79040581189389266572 

03.19 = -1.0153714219277431776 

03.20 = -0.25660116699213786317 

03.21 = 0.24235775601107244412 

03.22 = 0.032202698057872206652 

03.23 = -0.026566452015385539718 

04.1 = 0.0085760788432811750313 

04.2 = 0.013512845571180982502 

04.3 = 0.022567253340913112846 

04.4 = -0.058307759363885212745 
04,6 = -0.12041195448680664419 
04,6 = 0.22948717483369626712 



04,7 = 0.54863763-561981.341312 
a4,8 = 0.41924435198770363142 
04,9 = 0.25222427818457038144 



04,10 




-0.029868530356228792693 


04,11 




-0.30540634608037589943 


04,12 




-0. 18020289186124353787 


04,13 




-0.015353428793734303657 


04,14 




0. 10180220964078586702 


04,15 




0.23447509646113431046 


04, 16 




0.040929410205375688302 


04,17 




-0. 17564595949066362349 


04,18 




-0.051894216445801530355 


04,19 




0.065022578066260344054 


04,20 




0.015872184741392310993 


04,21 




-0.014923762934166461295 


04,22 




-0.0019211515134671124718 


04,23 




0.0015849038302656338915 
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The interior of If2C is given by: 
as = -0.0087234986549392520415 
ar = 0.010574245783011317441 
as = 0.041177278191813548941 
as = -0.037706537233780062173 
04 = -0.090162662747982636209 



03 = 0.018424137003272281874 
02 = 0.12963112574463663703 

oi = 0.25870815444749646286 
oo = 0.35615551493294340455 
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